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Abstract: The errors of low-cost inertial sensors, especially Micro-Electro Mechanical 
Systems (MEMS) ones, are highly dependent on environmental conditions such as the 
temperature. Thus, there is a need for the development of accurate and reliable thermal 
compensation models to reduce the impact of such thermal drift of the sensors. Since the 
conventional thermal calibration methods are typically time-consuming and costly, an 
efficient thermal calibration method to investigate the thermal drift of a full set of gyroscope 
and accelerometer errors (i.e., biases, scale factor errors and non-orthogonalities) over the 
entire temperature range in a few hours is proposed. The proposed method uses the idea of 
the Ramp method, which removes the time-consuming process of stabilizing the sensor 
temperature, and addresses its inherent problems with several improvements. We change 
the temperature linearly for a complete cycle and take a balanced strategy by making 
comprehensive use of the sensor measurements during both heating and cooling processes. 
Besides, an efficient 8-step rotate-and-static scheme is designed to further improve the 
calibration accuracy and efficiency. Real calibration tests showed that the proposed method 
is suitable for low-grade IMUs and for both lab and factory calibration due to its efficiency 
and sufficient accuracy. 

Keywords: MEMS inertial sensors; IMU; thermal calibration; turntable; temperature chamber 
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1. Introduction 

Integrated Global Navigation Satellite System (GNSS)/Inertial Navigation Systems (INS) systems 
are becoming widely used techniques to provide precise and reliable navigation information (i.e., 
position, velocity, attitude, etc.) in fields such as surveying and mapping, spatial information systems 
and intelligent transport [1]. A GNSS receiver is mainly used to provide high absolute position 
accuracy with long-term stability in ideal conditions, but has certain limitations in challenging areas 
(e.g., city downtowns, inside tunnels, under heavy tree canopies, etc.). On the contrary, as fully 
self-contained systems, INS measures the specific forces and angular rates by inertial sensors [i.e., 
accelerometers and gyroscopes (gyros)] and determines the motion of a body with respect to an inertial 
frame of reference [2]. Due to the integration of the inertial measurement, the sensor errors will 
accumulate and grow, resulting in increasing position and velocity errors. Since GNSS and INS are 
complementary systems, they can be integrated to give a system that has several advantages over 
each individual component separately. In such integration, the GNSS-derived positions and velocities 
are updating information through a Kalman Filter (KF) while the IMU is used to provide the 
navigation information during GNSS signal outages and for fast GNSS signal reacquisition [3]. 
However, traditional INS devices are bulky, expensive and complex, which often precludes their use 
outside specialized application fields, especially in civilian areas [4]. 

During these decades, advances in Micro-Electro-Mechanical Systems (MEMS) technology 
combined with the miniaturization of electronics have made it possible to produce chip-based inertial 
sensors. These MEMS chips have become ideal candidates for various applications since they are 
small, light-weight, low power and are extremely low-cost and reliable [5]. At the same time, the cost 
reduction of GNSS receivers has also promoted the development of low-cost navigation techniques. 
Accordingly, low-cost GNSS receiver-chip/MEMS INS integrated systems are becoming widely-used 
navigation techniques. 

However, the errors of the low-cost sensors, especially the MEMS ones, will change with time and 
are highly dependent on environmental conditions such as temperature [6]. Therefore, even though 
in-lab calibration at room temperature is known to be a useful way to remove the major part of the 
deterministic sensor errors, the actual values of the sensor errors can be different from those obtained 
through calibration processes due to the difference between the operational and calibration 
temperatures. In real tests, we find that the IMU sensor errors, not only biases, but also scale factors 
and non-orthogonalities, may vary significantly with temperature. These changes reach thousands of 
deg/h and many thousands of \ig for gyros and accelerometers biases, and thousands of ppm for the 
sensor scale factors and non-orthogonalities. Such errors, if not compensated for, will accumulate and 
lead to attitude and position drifts [7]. Approximately, for 2-D navigation, uncompensated gyros and 
accelerometers biases result in position errors of (l/6)b g (T)gt 3 and (l/2)b a (T)t 2 , respectively [5], 
where b g (T) and b a (T) are the gyro and accelerometer biases at temperature T\ t is the time that INS 
work alone; g is the local gravity value. Also, the uncompensated sensor scale factor errors may 
introduce position errors proportional to time squared. Thus it is important to do thermal calibration for 
MEMS IMUs so as to develop accurate and reliable thermal models. 

However, the conventional thermal calibration methods are typically time-consuming and costly. In 
this paper, an efficient thermal calibration method to investigate the thermal drift of a full set of 
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gyroscope and accelerometer errors (i.e., biases, scale factor errors and non-orthogonalities) over the 
entire temperature range in a few hours is proposed. 

2. Previous Works 

2.1. Calibration and Modeling of Sensor Errors 

The inertial sensor errors can be divided into two types: deterministic (systematic) errors and 
stochastic errors. Calibration is an important way to remove the major part of the deterministic IMU 
sensor errors. Various IMU calibration methods have been proposed according to different grades of 
IMUs and different applications [2,8-12]. On the contrary, the stochastic errors can only be modeled 
instead of being calibrated and removed. Among various stochastic modeling methods, Allan Variance 
is commonly used to determine the characteristic of the underlying random processes [13]. 

The thermal drifts of the sensor errors are normally regarded as deterministic errors. To reduce the 
thermal drift of the sensor errors, the following two processes are needed: (1) Thermal calibration: to 
develop accurate and reliable thermal models of the sensor errors, i.e., build the relationship between 
the sensor errors and the sensor temperature; (2) Thermal compensation: to compensate the thermal 
drift of the sensor errors according to their online temperature during the operation process of the 
IMUs. Both these processes are dependent on the temperature provided by the internal temperature 
sensors of the IMUs. 

Official thermal calibration needs professional equipment such as a thermal chamber and a turntable 
(optional). Such professional equipment can be used to get reliable thermal calibration results over a 
large temperature range. However, when using these high-cost equipments, it is also important to 
improve the calibration efficiency, so as to reduce the equipment cost. 

The purpose of thermal calibration is to determine the IMU sensor errors under different 
temperature points. Therefore, thermal calibration is commonly used based on the following 
assumption: the thermal drift of a sensor error is only related to the temperature of the sensor core. 
That is, corresponding to a certain temperature, the sensor errors will have the same values. Based on 
this assumption, the relationship between the sensor errors and the temperature can be built directly 
through either curve fit (as used in this paper) or a lookup table, if the corresponding sensor errors at 
different temperatures are known. 

There are currently two main approaches for thermal calibration: the Soak method and the Ramp 
method [2]. The Soak method works on the premise of stable sensor temperature while the Ramp 
method works based on time- varying sensor temperature. 

2.2. Temperature Soak Method 

The Soak method calibration works as follows: (a) Stabilize the temperature of the sensors or IMUs 
at a certain temperature point; (b) Record the sensor measurements and calculate the sensor errors 
corresponding to such point; (c) Repeat this process at several typical temperature points, then a series 
of sensor errors along with temperature data can be obtained. After that, the IMU sensor errors at other 
temperature points can be calculated by interpolation [14]. The Soak method is used for both the 
calibration of deterministic errors and the modeling of stochastic errors [15-21]. By stabilizing the 
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temperature, the Soak method can provide most reliable values of the sensor errors at the chosen 
temperature points. The main fact which restricts the use of the Soak method is the time required and 
the corresponding costs related to the equipment, time and manpower. The process of stabilizing the 
sensor temperature inside a MEMS IMU is usually time-costing. 

2.3. Temperature Ramp Method 

Different from the Soak method, the Ramp method calibration works as follows: (a) Control the 
temperature of the thermal chamber while it is continuously linearly increased or decreased to cover 
the full temperature range; (b) During this process the IMU measurements are recorded and the IMU 
sensor errors are calibrated [22,23]. The Ramp method is fundamentally faster in principle because it 
removes the stabilization process, and works based on linearly changed sensor temperature. However, 
there are two main inherent problems of the Ramp method, especially when it is used for the 
calibration of the IMUs, instead of the sensors: 

Issue #7: It is not realistic to get the real temperature of the inertial sensors: there are temperature 
differences between the inertial sensors and the temperature sensors inside the IMUs because they 
have different cores. What we can get is the temperature at the core of the temperature sensor, not that 
of the inertial sensors. Moreover, such temperature differences will be different under various 
temperature changing conditions. This issue is similar to the hysteresis effect. 

Issue #2: The chamber temperature changes during a calibration scheme (one set of IMU motions, 
with which the IMU outputs can be used to calculate for a set of IMU sensor errors; during the whole 
temperature calibration process, the turntable repeats the calibration scheme at different temperatures). 
This will lead to the changes of the IMU sensor errors to be calibrated, which in turn cause the 
calibration errors because the raw sensor data is not collected at the same temperature. 

Both these two problems can impair the performance of the Ramp method. This can also explain 
why the process of stabilizing the temperature is needed in the Soak method. In this paper, we use the 
idea of the Ramp method, which removes the time-consuming process of stabilizing the temperature. 
Also, several ways to ease its inherent problems are considered. The details of the proposed method 
are shown in the next section. 

3. Methodology 

This section includes the basic idea of the proposed thermal calibration method, including the 
ways to deal with the inherent issues of the Ramp method, the calibration equipment used and the 
corresponding calibration scheme. 

3.1. Basis Idea 

We use the idea of the Ramp method, but consider the following ways to mitigate its inherent 
problems for improvement: 

Improvement #7: take a balanced heating-and-cooling strategy to compensate the thermal 
calibration errors. To be specific, we change the temperature linearly for a complete cycle, and make 
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comprehensive use of the sensor measurements under both the linear heating and cooling processes at 
the same changing rate to get the final thermal calibration results at the same temperature point. 

Improvement #2: use the same calibration scheme during the whole calibration process. This can 
compensate the errors caused by the temperature changes during each calibration scheme to some 
extent after considering the heating and cooling processes comprehensively. In another word, the 
impact of the temperature changes in the calibration scheme during the heating process can be 
somehow cancelled with the impact during the cooling process at the same temperature point. 

Improvement #3: narrow down the temperature change during each calibration scheme. One 
possible way for this is to calibrate the IMU sensor errors within the shortest period while the 
temperature is changing. We designed an efficient 8-step rotate-and-static scheme, which can calibrate 
the gyro and accelerometer biases, scale factors and non-orthogonalities of MEMS IMUs within 
only 2.5 min. 

The corresponding explanations to these improvements can be stated as below: 

Explanation #7: under the heating and cooling processes at the same temperature changing rate, the 
effects of the corresponding temperature differences between the inertial sensors and the temperature 
sensors can cancel each other approximately (i.e., they have the similar value but opposite signs). 

Explanation #2: using the same calibration scheme means that the impact of the temperature change 
at each sensor error can be cancelled to the maximum extent. For example, assume it is under heating 
process and the nominal temperature of the calibration scheme (i.e., the average temperature during a 
scheme) is 7\; when the IMU is rotated around the x-axis gyro clockwise, the temperature is T\ + &T\. 
Then, for the cooling process at the same temperature changing rate at 7\, the IMU rotation around the 
x-axis gyro clockwise will happen at the temperature of T\ - &T\. Therefore the impacts of the 
temperature variations from the nominal temperature of the calibration scheme can be cancelled by 
using the heating and cooling results. 

Explanation #3: If the calibration scheme is well designed to be more efficient, one set of 
calibration motions can be finished in a shorter period and the temperature change during the scheme 
will be smaller. Therefore the disturbance of the temperature change to the calibration result at the 
temperature point can be mitigated accordingly. 

Both Improvement #1 and #2 are utilized to mitigate Issue #7, while Improvement #2 and #3 aim to 
ease Issue #2. With such improvements, the inherent problems of the Ramp method can be mitigated. 

3.2. Calibration Scheme 

A dual-axis position turntable equipped with a thermal chamber (see Figure 1) located at Wuhan 
University is used to calibrate the IMU sensor errors over a wide temperature range. The performance 
characteristics of both the turntable and the thermal chamber are shown in Table 1 . 

The turntable is accurate enough to provide reference for the calibration of low-grade IMUs. 
Besides, there are high-precision locating pins on the mounting plate of the turntable to keep the IMUs 
axes align with the turntable axes. 
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Figure 1. Thermal calibration equipment. 




Table 1. Performance of thermal calibration equipment. 



Turntable 



Principal axis rotation range 


Continuous infinite 


Principal axis angular position accuracy 


±5" 


Tilting axis rotation range 


±95° 


Tilting axis angular position accuracy 


±5" 


Non-orthogonalities between axes 


±5" 


Thermal Chamber 


Temperature range 


-55 - +100 °C 


Temperature change rate 


±0.1 ~+5 °C/min linear 



Considering both the characteristics of the equipment and the calibration efficiency, the following 
8-step calibration scheme is designed (see Figure 2). This scheme allows the IMUs to experience all 
necessary orientations and rotations with minimum calibration actions. 

In each calibration scheme, each sensing axis of every sensor is pointed alternately up and down 
precisely (except z-axis due to the mechanical structure of the turntable), and the IMUs are rotated 
around each gyro axis both clockwise and counter-clockwise at accurately known angles. There are in 
total eight static positions and eight rotations. 

In our test, each static position takes 10 s and the rate of rotation is about 10 deg/s, i.e., about 9 s for 
each rotation. Due to such settings, the calibration scheme is highly efficient, and can provide enough 
calibration accuracy for low-grade IMUs, as shown in Section 3.3 for the calculation of the calibration 
errors caused by sensor noises. Each calibration scheme takes only about (10 + 9) x 8 = 152 s to estimate 
a full set of IMU sensor errors. 
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During the whole calibration process, the thermal chamber keeps on changing its temperature 
linearly. At the same time, the turntable repeats the calibration scheme in Figure 2 again and again. 
Due to the efficiency of the scheme, we obtain abundant data of IMU sensor errors together with 
temperatures over the whole temperature cycle. 

Figure 2. IMU actions in an 8-step calibration scheme. 
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3.3. Sensor Error Models and Calibration Computation 

The error models for the accelerometers and the can be written as 

f=b fl+ S fl f + N fl f + v fl 



& = \> g +S g co + N g co + v g 



(1) 

(2) 



where f and g> are the error vectors of the accelerometer-derived specific forces and the 
gyro-derived angular velocities, f and co are the true specific forces and the true angular velocities, b a 
and b g are the biases of the accelerometers and the gyros, S a and S g are the diagonal matrices 
containing the scale factor errors, N« and N g are the skew- symmetric matrix containing the 
non-orthogonalities, v a and v g represent the a ccelerometer and gyro noises. 

In each calibration scheme, the IMU sensor errors, including biases, scale factor errors and 
non-orthogonalities, are determined. The IMU data is sufficiently used: the static data (i.e., the IMU 
outputs at eight static positions) is used to calibrate the gyro biases and all the accelerometer errors, 
while the dynamic data (i.e., the IMU outputs during eight rotation steps) is used to calibrate gyro scale 
factors and non-orthogonalities. The accelerometer errors can be estimated by the least-square 
method [2], while the gyro errors are estimated through a two-step method [24]. The details of 
calibration computations are shown in Appendix A. 
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The calibration errors caused by sensor noises can be calculated as follows [25]. Since there are 

ARW 

eight static positions, gyro noises can lead to calibration error of gyro biases <J b _ noise = —= — , (la), 

V8 •>/'.» 



static 



where ARW are the angular random walk of gyros, Static is the time for each static position. Similarly, 
the calibration error of accelerometer errors caused by accelerometer noises can be calculated. For the 
calibration of gyro scale factors and non-orthogonalities, there are at least two rotations (positive and 
negative) for each gyro. Thus the gyro noise can cause calibration error of gyro scale factors or 

non-orthogonalities a _ . = ^ ^ rotate - , where Rotate is the time for each rotation. For example, 

g yj 2 • 90P 

assuming the gyro ARW is 0.05 deg/ s I ^///z (refer to that of typical MEMS gyros), the errors caused 
by gyro noise are about 20 deg/h for gyro biases, and about 1,200 ppm for gyro scale factors or 
non-orthogonalities, which is much smaller than the thermal drifts of the sensor errors of low cost 
MEMS gyros. 

3.4. Establishment of the Thermal Variation Models 

Based on the calibrated IMU sensor errors at different temperatures, the curve fit method was used 
to determine the continuous global temperature compensation models (i.e., the relationship between 
the sensor errors, i.e., biases, scale factors and non-orthogonalities, and the temperature). A 3rd-order 
polynomial is used to fit the calibration results obtained from each scheme: 

W(T) = A 0 +A i .T + A 2 .T 2 +A 3 .T 3 ( 3) 

where W(T)is the calculated value of the sensor errors, T is the IMU core temperature in real 
operation process, Ao, Ai, A2 and A3 are the fitting parameters. The established thermal model can be 
used to calculate the thermal drifts of the IMU sensor errors at each temperature and then remove such 
thermal drifts from the IMU outputs. 

4. Thermal Calibration Tests and Results 

4.1. Temperature Change Strategy 

To evaluate the accuracy of the proposed method, we carried out the tests of the conventional Soak 
method to get the reference. The IMU calibration results at more than 10 temperature points were 
collected to provide accurate reference value, even though the Soak method does not need so many 
points and so much time for common uses. We also designed three different chamber temperature 
profiles to test the repeatability of the proposed method. Table 2 describes the detail of the chamber 
temperature changes. Figures 3 and 4 show the temperature profiles and how the internal temperatures 
of the two tested IMUs (Xsens MTi-G and NV-IMU100) changed. In both Figures 3 and 4, the dashed 
lines are the chamber temperatures, while the solid lines are the temperature provided by the internal 
temperature sensors of the two tested IMUs. It is clear to see the temperature lags of the IMU cores. 
Besides, Figure 3 shows that the IMU temperatures were higher than the chamber temperature after 
stabilization due to the heat generated inside the IMUs. 
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Table 2. Description of temperature changes in tests of both Soak and proposed method. 



Method 



Chamber Temperature Changes 



Soak method 

(as reference: provides the 
best reference of IMU 
sensor errors we can get) 



Proposed method 



Profile 1 
Profile 2 
Profile 3 



Consists of two parts: 

Part 1 (heat-and-stay part): increased from -30 °C to +70 °C using 10 °C steps 
(i.e., step at -30, -20, -10, 0, 10, 20, 30, 40, 50, 60 and 70 °C). Kept fixed at each 
step for 1 h. 

Part 2 (cool-and- stay part): decreased from +70 °C to -30 °C using 10 °C steps. 
Also 1 h per step. 

Changed linearly at +2 °C/min for a complete cycle over -35 °C ~ +75 °C. 
Kept the temperature fixed for a period at both -35 °C and +75 °C to expand the 
temperature range inside the IMUs to nearly -10 °C ~ +70 °C. 
A full calibration cycle took about hours. 

Used the following three profiles to test the repeatability of the proposed method. 
25 °C — -35 °C (keep 30 min) — +75 °C (keep 60 min) —35 °C (keep 30 min) 
25 °C — +75 °C (keep 30 min) — -35 °C (keep 30 min)^+75 °C (keep 30 min) 
25 °C ^+75 °C (keep 40 min) — 25 °C (keep 20 min) —35 °C (keep 30 min) 
— 25 °C (keep 30 min) 
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Figure 3. Temperature profile in Soak method test. 
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Figure 4. Temperature profiles in proposed method tests. 
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4.2. IMUs Used in the Tests 

Two MEMS IMUs, i.e., Xsens MTi-G [26] and NV-IMU100 [27], were chosen to be the 
representative of low-grade IMUs to be tested. Their photos in the test and their characteristics are 
illustrated in Figure 5 and Table 3, respectively. Here we used the raw sensor data of MTi-G (without 
factory compensation). The NV-IMU100 has three vibrating ring gyros, and its output data had been 
partly compensated after factory calibration. 

Figure 5. Two tested MEMS IMUs and the installation on the mounting plate of 
the turntable. 




Table 3. Characteristics of tested IMUs. 



Characteristics 


MTi-G 


NV-IMU100 


Data Rate 


100 Hz 


166 Hz 


Operating Range 


-40-85 °C 


-40-80 °C 


Gyro Bias (lo ) 


3,600 deg/h 


4,000 deg/h 


Gyro White Noise (ARW) 


3.0 deg/Vh 


1.8 deg/Vh 


Gyro Scale factor error (lo ) 




5,000 


Accel. Bias (lo ) 


2,000 jig 


9,000 \ig 


Accel. White Noise (VRW) 


0.002 m/s2/VHz 


0.01 m/s2/VHz 


Accel. Scale factor error (lo ) 


3,000 ppm 


2,000 ppm 



4.3. Calibration Results of MTi-G 

Figure 6 shows the results of the biases and scale factor errors of gyros with the Soak method, as 
the example results of this conventional method. 

The thin solid lines in Figure 6 show that even when using the Soak method, there were some 
non-repeatabilities between the results during Part 1 (from -30 to 70 °C with 10 °C step) and Part 2 
(from 70 to -30 °C with 10 °C step). After checking the temperature data of the IMUs, we found that it 
is because the sensor temperature inside the IMUs had not been fully stabilized even after the chamber 
temperature outside the IMUs had been stabilized for 1 h. This means that more than 1 h is needed for 
each temperature point to get the ideal result of the Soak calibration method for MTi-G. However, 
we cannot afford more than 1 h (e.g., 2 h) for each point for in-lab tests. Thus, we finally used 1 h per 
point, but fitted the results of Part 1 and Part 2 (i.e., heating and cooling profiles) at the same 
temperature point to calculate the final IMU sensor errors at this temperature. This is the best result we 
can get and was regarded as the reference for the evaluation of the proposed method. 
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Figure 6. Thermal calibration results of MTi-G gyro biases and scale factor errors with 
Soak method. The thin solid lines are the results under separated heating and cooling 
processes (rightward and leftwards arrows indicate heating and cooling processes, 
respectively). The thick lines are the final curve fitted result (3rd-order polynomial fit). 
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Figure 7 shows the results of gyro biases and scale factor errors with the proposed method, as the 
example results of this new method. 

The raw results of the proposed method (the thin dashed lines) have some non-repeatabilities during 
heating and cooling processes. This might be caused by different temperature change rates. However, 
after fitting the raw results, we got the final result (the thick dashed line), which fit well with the 
reference result in Figure 6. 



Figure 7. Thermal calibration results of MTi-G gyro biases and scale factor errors with 
proposed method. The thin dashed lines show the raw calibration results (i.e., the result of 
each calibration scheme). Rightward and leftwards arrows indicate heating and cooling 
processes, respectively. The thick dashed lines are the final fitted results (3rd-order 
polynomial fit). 
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To test the repeatability of the proposed method, we run the chamber temperature as the 
three different profiles described in Table 2. Figure 8 shows the final results of the three tests of the 
proposed method with different profiles (three dashed lines) as well as that of the Soak method 
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(the solid line). Figure 8a-f show the curves of biases, scale factor errors, and non-orthogonalities of 
gyros and those of accelerometers. 

Figure 8. Thermal calibration results of MTi-G. In each subplot, three dashed lines 
correspond to the results of the proposed method with three different temperature profiles, 
and the solid line is the result of the Soak method. 
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It is clear that the sensor errors of MTi-G varied significantly with temperature, especially the gyro 
biases and scale factors. Over the entire temperature range, the changes reached 2,000 deg/h, 
10,000 ppm for MTi-G gyro biases and scale factor errors, and 3,000 \±g and 400 ppm for 
accelerometer biases and scale factor errors. There were also changes of about 5,000 ppm and 
1,000 ppm for gyro and accelerometer non-orthogonalities, which were beyond our intuition. These 
curves indicate that thermal calibration is important not only for biases and scale factors, but also for 
non-orthogonalities, at least in the case of MEMS IMUs. 

To evaluate how much the final results provided by the proposed method (the three dashed lines) fit 
with that of the Soak method (the solid line), we picked some typical temperature points and calculated 
the Root Mean Square (RMS) errors (using the results of the Soak method as reference) of the 
three temperature profiles, as shown in Table 4. The first 12 rows are the results of the gyro biases, 
scale factor errors and non-orthogonalities, while the last 12 rows are those of the accelerometers of 
MTi-G. Columns 2-6 show the RMS of the sensor errors under typical temperature points and column 
7 picks up the max values among them. 

Table 4. Errors (RMS) of MTi-G thermal calibration at typical temperature points 
(using the results of the Soak method as reference). 
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Over the entire temperature range, their RMS with respect to the result of the Soak method reached 
360 deg/h, 1,500 ppm and 1,300 ppm for gyro biases, scale factors and non-orthogonalities; and 
1,200 |ng, 100 ppm, and 400 ppm for accelerometer biases, scale factor errors and non-orthogonalities. 
The accuracy of the proposed method is shown in Table 5. Also, the original thermal drift levels of 
each sensor error over -10 ~ +70 °C are shown to make a comparison. 

The overall accuracy levels are much lower than the levels of the original thermal drifts of the 
MEMS sensor errors. The proposed calibration method has enough accuracy to significantly reduce the 
original temperature variations of the MEMS IMU sensor errors. 

The thermal calibration results of MTi-G also indicate that not all the sensors errors can be reduced 
significantly through calibration, since the errors of each sensor are randomly distributed in the range 
of its specifications and some errors may be small originally. What we can expect is that through 
calibration and compensation, the overall levels of IMU sensor errors can be reduce to a much 
lower level. 



Table 5. Calibration accuracy of proposed method (reference to Soak method) and the 
original thermal drift levels of the IMU errors over the whole temperature range 
(with MTi-G). 



Sensor Errors 


Accuracy of Thermal Calibration 


Original Drifts over-10 - +70 °C 


Gyro biases 


360 deg/h 




2,100 deg/h 


Gyro scale factor errors 


1,500 ppm 




11,000 ppm 


Gyro non-orthogonalities 


1,300 ppm 




6,000 ppm 


Accel, biases 


1,200 jig 




3,800 fig 


Accel, scale factor errors 


100 ppm 




400 ppm 


Accel, non-orthogonalities 


400 ppm 




1,000 ppm 



4.4. Calibration Results ofNV-IMUlOO 



Similarly, we calibrated the thermal drift of NV-IMU100 with both the proposed method 
(using three different profiles) and the Soak method. Figure 9 shows the results provided by both the 
proposed method (the three dashed lines) and the Soak method (the solid line). 

There were small gyro biases (less than 15 deg/h) due to the symmetry of vibrating ring gyros 
and the factory compensation. However, the change for gyro scale factor (e.g., y-axis) exceeded 
10,000 ppm. These characteristics matched the inherent feature of the vibrating ring gyros. 

Like Table 4, Table 6 shows the RMS errors (using the results of the Soak method as reference) of 
the proposed method at some typical temperature points. 

Over the entire temperature range, the RMS errors reached 5 deg/h, 2,200 ppm and 900 ppm for 
gyro biases, scale factor errors and non-orthogonalities; and 4,200 |ug, 700 ppm, and 600 ppm for 
accelerometer biases, scale factor errors and non-orthogonalities, respectively. The accuracy of the 
proposed method is shown in Table 7 as well as the original thermal drift ranges of each sensor error 
over -10 ~ +70 °C. The overall accuracy levels are much lower than the levels of the possible thermal 
drift of the MEMS sensor errors. 
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Figure 9. Thermal calibration results of NV-IMU100. In each subfigure, three dashed lines 
correspond to the results of the proposed method with three different temperature profiles, 
and the solid line is the result of the Soak method. 
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Table 6. Errors (RMS) of NV-IMU100 thermal calibration at typical temperature points 
(using the results of the Soak method as reference). 



Cnri cot* rrnrc 


—10 T 


10 T 


^0 T 


50 T 


70 T 


IVfdY 


D_gx (aeg/nj 


c 

J 


1 
1 


1 
1 


A 
U 


9 
z 


C 

J 


o_gy i^aeg/n; 


u 


1 
1 


1 
1 


9 
z 


a 


a 


D_gz ^aeg/n; 


a 


9 


1 
1 


A 
4 


A 
4 


A 
4 


os gx (ppm; 




JOJ 


^98 
JZo 


79 

/Z 




JOJ 


os _gy ippmj 


Z, 1HO 


z,ujy 


1 , D 1U 


^79 

J) / Z 


1 on 


9 1A^ 
Z, 14J 


5s_gz (ppm) 


/ 


/I ^7 
40 / 


J 44 




7Q9 

lyZ 


7Q9 

lyZ 


n_xy (ppm) 


Jyo 


^71 

J / 1 


AIL 
4 / j 


A^n 

4JU 


^78 

J /o 




n_xz (^ppmj 


A(\ 


09 

yL 


44 


1 A1 
141 


^7 
J / 


1 A1 
141 


n_yx ^ppmj 




^99 

JZZ 


/LQ7 
4^ / 


JJO 


AH^ 
4UJ 


27A 
0 /4 


n_yz (ppmj 


J / 


41 


9H 
ZU 


A7 
4 / 


1 1 J 


1 1 J 


n_zx (ppm) 


1UJ 


ol 


/in 

4U 


A9 
4Z 




1UJ 


n_zy (ppmj 




J 1 


4U 


A9 
4Z 


JU 


1UJ 


k\ 0 a/" / 1 1 r^r 1 

D_ax VJLXgJ 


4, 1 J D 


'X ^9H 


^ A1 7 
J, 41 / 




1 ^^7 


4, 1 j j 


r\ 0 "\ t 1 1 1 rr 1 

D_ay ^(igj 


1 ,OJ4 


1 ,jOO 


7HA 
/U4 


Q9A 

^/Z4 


9 ^9^ 

Z, JZJ 


9 ^9^ 

Z, JZJ 


D_az ^|ig; 




9 1 7H 

Z, 1 /U 


1 ^1A 
1,J 14 


1 , 1U1 


j,oOU 


J,oOU 


An O AT" / VMAtTl 1 

os ax yjpmj 




^9 
OZ 


AO 
4^ 


Do 


on 


yj 


8 s ay (ppm) 


1 1 J 


^H7 


4jU 


1 9H 
1ZU 


9^^ 

ZOj 


ou / 


8s_az (ppm) 


179 


146 


135 


153 


226 


226 


m_xy (ppm) 


380 


443 


490 


523 


511 


523 


m_xz (ppm) 


223 


106 


41 


16 


16 


223 


m_yx (ppm) 


547 


398 


472 


521 


313 


521 


m_yz (ppm) 


169 


149 


124 


50 


134 


169 


m_zx (ppm) 


161 


83 


102 


87 


135 


161 


m_zy (ppm) 


205 


112 


100 


110 


69 


205 



Table 7. Calibration accuracy of proposed method (reference to Soak method) and the 
original thermal drift ranges of the sensor errors over the whole temperature range 
(with NV-IMU100). 

Sensor Errors Accuracy of Thermal Calibration Original Drifts over -10 - 
+70 °C 

Gyro biases 5 deg/h * 13 deg/h 

Gyro scale factor errors 2,200 ppm 10,000 ppm 

Gyro non-orthogonalities 900 ppm 1 ,000 ppm 

Accel, biases 4,200 jig 14,000 jig 

Accel, scale factor errors 700 ppm 1,800 ppm 

Accel, non-orthogonalities 600 ppm 600 ppm 

* The NV-IMU100 has three vibrating ring gyros, and the thermal drifts of its gyro biases had been partly 
compensated by factory calibration. 
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The proposed thermal calibration method can provide results close to those of the conventional 
Soak method, and can significantly reduce the calibration time. Thus it is especially useful for MEMS 
IMUs, which require low cost and efficient calibration. Please note that the thermal calibration 
accuracy not only depend on the calibration method itself, but also is affected by the stochastic errors 
of the calibrated inertial sensors. This can explain why the same method got different calibration 
accuracies for the two tested IMUs. For example, the calibration accuracy of gyro biases of 
NV-IMU100 can reach 5 deg/h because the original gyro biases are relatively stable. To further 
investigate the effect of thermal calibration, another kind of tests called thermal compensation tests 
were performed, as described in Section 5. 

5. Thermal Compensation Tests and Results 

5.7. Test Description 

The thermal compensation tests were designed to check how much the performance of the inertial 
sensors can be improved by compensating their outputs with the established thermal models. There 
were totally six tests, including three dynamic tests and three static tests, which are described in Table 8. 
For each test, the thermal chamber was utilized to provide different temperature environment (from 
-10 °C to +70 °C), and the turntable was used to provide precisely known reference. Every epoch of 
IMU outputs was corrected using the full set of IMU sensor errors provided by the thermal model at 
the same temperature (i.e., the thermal compensation). Then the effect of the thermal compensation 
was evaluated by comparing the output errors (i.e., the difference between the sensor outputs and the 
reference input) before and after the thermal compensation. To evaluate the performance of all the 
sensors before and after compensation, every accelerometer axis was kept on pointing downwards 
precisely for about an hour respectively. Besides, the IMUs were kept on rotating around each gyro 
axis at an angular rate of 10 deg/s, which was set according to the dynamic of land vehicles, for about 
an hour. The whole process took about six hours. In each test, we focused on the output error of one 
sensor. For example, we investigated the x-axis gyro output with a reference input of 10 deg/s in the 
dynamic test 1; and investigated the x-axis accelerometer output with a reference input of the local 
gravity in the static test 1 . The IMU outputs were averaged during each second to reduce the impact of 
the noises. The thermal models established through thermal calibration test 2 using the 2nd 
temperature profile were used in this test, since the results provided by test 2 seemed to have the 
largest differences comparing with those of the Soak method. The parameters of the 3rd order 
polynomial thermal models of applied in the compensation tests, i.e., A 0 , A\, A 2 and A 3 in Equation (3), 
for the sensor errors of the two tested IMUs are shown in Table 9. The test results of MTi-G and 
NV-IMU100 are shown in Section 5.2 and 5.3, respectively. 
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Table 8. Descriptions of the six thermal compensation tests. 



Motions of the IMUs (Provided by the Turntable) 



Dynamic test 1 


IMUs kept rotating around their x-axis at 10 deg/s. 


Dynamic test 2 


IMUs kept rotating around their y-axis at 10 deg/s. 


Dynamic test 3 


IMUs kept rotating around their z-axis at 10 deg/s. 


Static test 1 


IMUs kept static with x-axis pointed upwards. 


Static test 2 


IMUs kept static with y-axis pointed upwards. 


Static test 3 


IMUs kept static with z-axis pointed upwards. 



Chamber Temperature 



Before each test, the temperature was stabilized at -10 °C; 
During each test, the temperature was changed from -10 °C to 70 °C in one hour. 

Table 9. The parameters of the 3rd order polynomial thermal models of the IMU errors 
applied in the thermal compensation tests, i.e., Ao, Ai, A2 and A3 in Equation (3). The 
thermal models can be used to calculate the thermal drifts of the IMU sensor errors at each 
temperature and then remove such thermal drifts from the IMU outputs. 
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5.2. Compensation Results ofMTi-G 

Figures 10 and 11 show the compensation results of MTi-G. The blue and red lines are the output 
errors without and with the compensation using the established thermal models. 

It is clear that the original output error of the z-axis gyro was significantly larger and changed with 
some long-term trends. After compensation, such error was significantly reduced, and the long-term 
trends were eliminated. However, there was no significant improvement for x-axis and y-axis gyros, 
since their original output errors were relatively small. As for the accelerometers, the improvement 
was significant along x-axis and z-axis. However, the y-axis accelerometer got an output error of about 
2,000 |ng after compensation. As has been explained in Section 4.3, not all the sensors errors can be 
reduced significantly through the thermal compensation, because the thermal drift level of each sensor 
are randomly distributed in the range of its specifications and some may be small originally. The 
thermal compensation can reduce the overall thermal drifts of the IMU outputs to a much lower level. 



Figure 10. MTi-G gyro output errors with and without compensation. 
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Figure 11. MTi-G accelerometer output errors with and without compensation. 
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The statistical values (i.e., the mean and the RMS values) of the IMU output errors over the entire 
temperature range were calculated and shown in Table 10. Without compensation, the RMS values 
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were 279, 197 and 2,167 deg/h for gyros and 3,697, 1,097 and 11,130 |ug for accelerometers. After the 
compensation, the RMS values became 189, 254 and 418 deg/h for gyros and 2,342, 2,259 and 1,055 
|ng for accelerometers. Both the gyro and accelerometer along the z-axis and x-axis were improved 
significantly, while those along the y-axis became worse. Generally speaking, with the thermal 
compensation, the output errors of MTi-G were reduced to within 420 deg/h and 2,400 |ug for gyros 
and accelerometers over the entire temperature range. To make it clear, the RMS values were shown in 
Figure 12. The blue and red bars were the values without and with thermal compensation, respectively. 

Table 10. Statistical results of MTi-G output errors with and without compensation. 
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Figure 12. RMS values of MTi-G output errors with and without compensation. 
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5.3. Compensation Results ofNV-IMUlOO 

Similarly, Figures 13 and 14 show the NV-IMU100 output errors without and with 
thermal compensation. The NV-IMU100 has small gyro biases due to the symmetry of the vibrating 
ring gyros and the factory compensation. However, there were also gyro output errors caused by the 
scale factor errors and non-orthogonalities, especially for the y-axis gyro, whose scale factor varied for 
thousands of ppm over the entire temperature range. Consequently, the gyro output error reached 
several hundred deg/h when the input angular velocity was 10 deg/h. Clearly, such output errors were 
reduced significantly after the thermal compensation. The errors of accelerometers, especially the 
x-axis and y-axis accelerometers, were also reduced dramatically. The systematic errors in the outputs 
were corrected, and long-term trends in the outputs were eliminated. The statistical values of the 
NV-IMU100 output errors were shown in Table 11, and the RMS values were illustrated in Figure 15. 
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Figure 13. NV-IMU100 gyro output errors with and without compensation. 
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Figure 14. NV-IMU100 accelerometer output errors with and without compensation. 
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Table 11. Statistical results of NV-IMU100 output errors with and without compensation. 
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Figure 15. RMS values of NV-IMU100 output errors with and without compensation. 
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6. Conclusions 

An efficient thermal calibration method is proposed to investigate the thermal drift of a full set of 
IMU sensor errors (i.e., biases, scale factor errors and non-orthogonalities) over the entire temperature 
range in 4 h. Results of thermal calibration tests with two MEMS IMUs over the temperature range of 
-10 ~ +70 °C showed that the proposed method provided results close to those of the conventional 
Soak method, and significantly reduced the calibration time. The accuracy of the proposed method 
over the entire temperature range (reference to the Soak method) was shown in Table 12. 

Results of the separate thermal compensation tests showed that the performance of the inertial 
sensors can be significantly improved by compensating their outputs using the established thermal 
models. The max RMS values of the IMU output errors without and with thermal compensation were 
shown in Table 13. 



Table 12. Accuracy of proposed method (reference to Soak method). 



Sensor Errors 


Xsens MTi-G 


NV-IMU100 


Gyro biases 


360 deg/h 


5 deg/h * 


Gyro scale factor errors 


1,500 ppm 


2,200 ppm 


Gyro non-orthogonalities 


1,300 ppm 


900 ppm 


Accel, biases 


1,200 jig 


4,200 fig 


Accel, scale factor errors 


100 ppm 


700 ppm 


Accel, non-orthogonalities 


400 ppm 


600 ppm 



* The NV-IMU100 has three vibrating ring gyros with stable bias inherently, and the thermal drifts of its 
gyro biases had been partly compensated by factory calibration. 



Table 13. RMS of IMU output errors with and without thermal compensation. 



Sensors 


Xsens MTi-G 


NV-IMU100 


Compensated Uncompensated 


Compensated Uncompensated 


Max Gyro output error (deg/h) 
Max Accel, output error (ug) 


418 2,164 
2,342 3,697 


48 314 
5,410 17,356 



Due to the high efficiency and sufficient accuracy of the proposed thermal calibration method, it is 
especially useful for both factory and lab calibrations of low-cost IMUs (e.g., MEMS), which require 
low-cost and efficient calibration. It can also promote better utilizations of low-cost MEMS-based 
sensors and IMUs. 
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Appendix A. Calibration computation of the proposed calibration scheme 

Estimation of Accelerometer Errors: 

To estimate a full set of accelerometer errors, the output of a triad of accelerometers is represented 
in matrix form: 

'/; 

/, 

/, 
i 
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ax 


/, 




m xy 


S y 


m zy 
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ay 


I 






m yz 


S z 


b 

az 



M 



(a.l) 



The diagonal s elements are the scale factors, the off diagonal m elements represent the 
non-orthogonalities and the b components are the biases. Assume SI, S2, S3, S4, and S5 are the IMU 
positions with x-axis upwards, x-axis downwards, y-axis upwards, y-axis downwards and z-axis 
downwards. Since the IMU is perfectly aligned with the turntable axes through the location pins on the 
mounting plate, the given specific forces of the five positions can be represented as follows: 
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(a.2) 



Then the design matrix can be denoted by A and the measured specific force of the accelerometer is 
denoted by U: 



A = 



f f f f f 

l l l 2 l 3 l 4 l 5 



11111 

U = K u 2 u 3 u 4 u 5 ] 

In this case the column vector of the U matrix should be: 



(a.3) 
(a.4) 



I 



I 



(a.5) 



X -axis pointing down 



- -lx- axis pointing up 

U3, U4 and U5 are similar with ui and U2 Then the M matrix can be estimated by the 
least- square method: 

M = U- A r (A- A r ) _1 (a .6) 

Estimation of Gyro Biases: 

The gyro biases are calculated using static gyro outputs. 

Let CO; = \jb ix cb iy a> iz ~J i = 1,2,3,4,5 are the mean of gyro output angular velocities under static 
positions SI, S2, S3, S4, and S5. The elements of the vector represent the components along three axes. 
By averaging a couple of measurements, e.g., c& x and co 2 or G) 3 and o> 4 , the impact of the rotational 
angular velocity of the earth can be cancelled. The gyro biases are calculated as: 
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1 4 1 4 

b = d)j b = — / ( 



^ £=1 ^ fc=l 

When the z-axis of the IMU is pointed downward, the z-axis gyro bias can be calculated using the 
gyro output minus the reference input when the gyro points vertically (which is not depend on the 
north finding of the turntable): 

b gz =& Sx + a) e smq> ( a . 8 ) 
where co e is the rotational angular velocity of the earth, <p is the latitude. 

Estimation of Gyro Scale Factors and Non-orthogonalities: 

Let L. = ^L ix L iy L. z J / = 1,2,3,4,5,6 be the integration of the gyro outputs when the IMU is 
rotated counterclockwise and clockwise around x-axis, counterclockwise and clockwise around y-axis, 
and counterclockwise and clockwise around z-axis for the angle of L re f, respectively. There are two 
factors that eliminate the impact of the rotational angular velocity of the earth: (a) In the designed 
calibration scheme, when the IMU is rotated around one axis both counterclockwise and clockwise, 
this axis is pointed toward its specific direction (as shown in Figure 2). (b) The rotation data 
(counterclockwise and clockwise) with the same time length is used for the calibration computation. 

The gyro scale factors can be calculated as: 

s = kxZhx _ L 3y~ L 4y _ h Z ~h Z 

8X ^ ref ' s sy ~ ' ^ " (a.9) 

ref ref 

The value of L re f is 90 deg in our calibration scheme. The non-orthogonalities cause each axis to be 
affected by the signal of the other two axes. When the IMU is rotated around one axis, the gyro output 
of another axis will be affected by this rotation due to the non-orthogonalities between these two axes. 
The non-orthogonalities can be calculated as: 



Xy rs j , XZ rs j , }>X rs j , yZ rs j , ZX rs j , Zy rs j 

^^ref ^^ref ^^ref ^^ref ^^ref ^ref 

where ny is the non-orthogonalities of i-axis to j-axis. 



(a. 10) 
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